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Abstract 

We investigate gauge invariant cosmological perturbations in a spatially flat 
Friedman-Robertson- Walker universe with scalar fields. It is well known that 
the evolution equation for the gauge invariant quantities has exact solutions 
in the long- wavelength limit. We find that these gauge invariant solutions 
can be obtained by differentiating the background solution with respect to 
parameters contained in the background system. This method is very useful 
when we analyze the long-wavelength behavior of cosmological perturbation 
with multiple scalar fields. 
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I. INTRODUCTION 



The theory of gauge invariant cosmological perturbations [1-3] is an important tools 
when we investigate the origin and the evolution of the structure in our universe. In the 
context of general relativity, the treatment of linear perturbation of any matter field requires 
great care because these perturbations are in general not gauge invariant. We must evaluate 
gauge invariant quantities to extract physical meaning. 

One of the main application of the gauge invariant perturbation is the theory of infla- 
tionary universe, in which the gauge invariant treatment of perturbations is crucial. The 
inflation predicts that the seed of density fluctuation, which evolves to the structures in our 
universe, is generated as quantum fluctuation of the inflaton field. The created fluctuations 
is stretched by deSitter expansion and their wavelength exceeds the Hubble horizon scale. 
As the gauge dependence of perturbations becomes conspicuous on super horizon scale, we 
must use gauge invariant perturbation to handle such a situation. But the derivation of the 
gauge invariant equation and evaluation of its solution needs tedious calculation in general, 
especially for the model with multiple scalar fields which is considered in the context of 
hybrid or extended inflation and scalar-tensor theories. 

We here pay attention to the Einstein gravity with minimally coupled scalar fields (f>. 
As a background space, we assume a flat Friedman-Robertson- Walker(FRW) universe of 
which dynamical variable is a scale factor a = e a . We use a following gauge invariant 
combination(Mukhanov's variable) [4]: 



where 5(f) is a perturbation of the scalar field, y? is a perturbation of three curvature and 
H = a is a Hubble parameter. The evolution equation for Q is given by 



It is known that this equation has exact solutions in the long- wavelength limit k — > [5]: 



Q = 8(f) - —<p, 



(1) 
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ft H 2 



(3) 



HJ e 3a <j) 2 



where c±,C2 are arbitrary constants. Therefore if we can obtain the background solution 
(a(t) , (j>(t)) , it is possible to predict the behavior of the long-wavelength gauge invariant 
perturbation without solving the evolution equation of the perturbation. 

Our question is why we can have the exact solution for k — > (0-mode) gauge invariant 
perturbation. The 0-mode perturbation is a perturbation of homogeneous mode, it must be 
contained in a homogeneous background system. So we expect that the evolution equation 
for the 0-mode perturbation can be derived by the analysis of the mini-super space model 
which has no inhomogeneous mode. 

In this paper, we aim at clarifying the relation between the gauge invariant 0-mode 
solution and the background solution. We will show that the evolution equation and the 
solution of the gauge invariant 0-mode perturbation can be obtained within the perturba- 
tion of mini-super space model. In Sec. II, we use one dimensional autonomous system to 
demonstrate that the perturbed equation and its solution is obtained by differentiating the 
background solution with respect to parameters contained in the background solution. In 
Sec. Ill, we treat a mini-super space model with scalar fields and derive the 0-mode gauge 
invariant perturbation equation(Mukhanov's equation) and its solutions. We use Hamilton- 
Jacobi method which is used by several authors [6,7] to derive the evolution equation of the 
gauge invariant cosmological perturbation. This method does not need to specify any gauge 
condition during calculation, it is suitable to apply to gauge invariant perturbation. Sec. IV 
is devoted to summary. We use the unit in which c = H = 8nG = 1 and denote the partial 
derivative of a function F with respect to a some variable A by Fa throughout the paper. 



As a demonstration, we consider a one dimensional autonomous system with a Hamilto- 
nian 



II. TOY MODEL EXAMPLE 
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(4) 



Equation of motion is 



x + V'(x) = 0. 



(5) 



The solution of this equation is written as x = x(t + to) where to is an origin of time and an 
arbitrary constant. The equation for the perturbation is obtained by splitting x = x^ + 5x 
and linearize with respect to 5x: 



This equation has the same form as the equation of x which is obtained by differentiating 
the background equation (5) with respect to time variable t: (x)" + V"(x)x = 0. As x = x to , 
a solution of the perturbation equation is obtained by differentiating a background solution 
with respect to the parameter t . 

We explain this result more rigorously using Hamilton- Jacobi(H- J) method which is 
suitable to treat a cosmological model in next section. The H-J equation for a generating 
function S is 




(6) 




(7) 



and the evolution equation is 



(8) 



Eq.(7) is solved by S 



Et + W(x)(E is a separation constant): 



W 2 X + 2V(x) = 2E, 



x = W x . 



(9) 



The solution is 




(10) 
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where t is an arbitrary constant. From this we get x = x(t + t ; E) or E = E(x,t + t ). The 
background solution is specified by the two parameters (t , E) and the perturbed solution 
is defined by the small change of the constants (to + St , E + 5E). 

The generating function is written as S = + F where is a background part and 
F is a perturbed part which consists of quadratic terms of St and 5E. F satisfies 

^ + S^F x + ^ = 0. (11) 

Now we change variables from (x, t) to (t(x; E),t , E(x, t + to)): 

d dr d dt d dE d 
di = dtlh + ~dtdfo + ~dtdE 

= -k +E 'w < 12 > 

d dr d dt d dE d 
— = 1 ^ 1 

dx dx dr dx dt dx dE 

The Eq.(ll) becomes 

OF 1 / [OF dF\ „ dF\ 2 . . 

a7 + 2{ T '{-37 + wJ +E 'aE) = ' < 14) 

As we are interested in linear perturbation due to the change of the parameters (Sto,SE), 
we replace ^ — > tJ^, — > Using a new variable 5i = 5t + t e 5E and taking into 
account the quadratic form of F, the H-J equation becomes 

dF 1 ( dF\ 2 , . 

* + 2i5(«j = »• (15) 

From this equation, we get a Hamiltonian for the perturbation variable 5t: 

^ = ^> (16) 
where P$ t is a conjugate momentum of St. The equation of motion is 

St + 2-St = 0. (17) 

x 

The equation for a variable Sx = = £ fo <5t + x E SE becomes 



5x + V"(x)5x = 0. (18) 

This is our desired result The two independent solutions are given by 

5x = x to ,x E - (19) 

The explicit form of the solution is obtained by differentiating Eq.(10) with respect t 
and E. Regarding x as a function of (t + t , E) and differentiating Eq.(10) with respect to 
t and E, 

d_ J_ 



and we have 

xt = W X = x, 

TTr f x dx . r* dt 

"-J n; ' / (21) 

The perturbed solution is obtained by differentiating the background solution with respect 
to the parameters which specify the background solution. 

III. GAUGE INVARIANT LONG- WAVELENGTH SOLUTION FROM 
MINI-SUPER SPACE SOLUTION 

A. a single scalar field case 

We consider a spatially flat FRW model with a scalar field. Dynamical variables are a 
scale factor a = e a and a scalar field <fi. The Hamiltonian for this system is 



H = Ne 



-3a 



-^Pl + \pl + e* a V{<P) 



(22) 



where iV is a lapse function, P a and are conjugate momentum of a and <fi, V(<f>) is an 
arbitrary potential for the scalar field. The equations of motion are 
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a 



Ps = 



on 

dK = 
_&H 

da 
&H 



N 
~6 



_l_ p -3ap 



dU 

dK 



Ne 



-3a 



Ne~ 3a 



(23) 



d(j) 



= Ne ia V'((j)). 



The Hamilton- Jacobi equation follows from the Hamiltonian constraint with the canonical 
conjugate momentum replaced by derivative of a generating function S, P a 



dS p 



as. 

d<t>- 



"§ + f + e 6a V(<j>)=0. 



(24) 



The evolution equations are 



a 
N 



,-3a 
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N ~ V 

The background solution is obtained by assuming = — 2e 3Q if (0). Then 

- 3H 2 + 2H\ + V((j>) = 0, 

- = H, 
N 



(25) 
(26) 



N 



-2m. 



(27) 
(28) 

(29) 



The solution of Eq.(27) has a constant of integration C and the solution can be written 
H = H((fi; C). By differentiating Eq.(27) with respect to C and using the evolution equation 
(28) and (29), we have 



H c = De- 3a = e - 3(a+ao) , 



(30) 



where D and a are arbitrary constants. From this, we can express a + a as a function of 
(0,C): 



a + «o = <5(0; C). 



(31) 



If we invert this relation, we get = 0(a + a$, C). «o is the origin of the scale factor and C 
determines the initial value of the scalar field 0. The background solution is characterized 



by two parameters (ceo, C) and perturbation is defined by small change of these parameters 
(«o + C + 5C). 

We split the generating function into the background and the perturbed part: S = 
S<® + F. F satisfies 

C(0) p2 p2 

--^-^ + ^- I f + f = 0. (32) 

Change variables from (a, 0) to (a, cto, C): 

d da d da d dC d 
da da da da da da dC 

° ' gJL (33) 



da dC 
d da d da® d dC d 
d<fi d(f) da d(f) da d<p dC 

~ ( 9 9 \ n 9 

= a *{m + da- ) +c *dc- (34) 

Then the H-J equation (32) becomes 

3aTT dF 1 / dF „ dF\ 2 1 / fdF dF\ „ dF\ 2 n 
G H lm-T2{-d^ + C °dc) + 2{^{da + da- ) +C ^dc) =°- (35) 

To extract the linear perturbation part, we replace ^ — >• q^,^ —> -^q- Using a new 
variable 1Z = 5ao + acSC and assuming a quadratic form of the generating function F with 
respect to 7Z, we have 

dF H (dF\ 2 n 

Ba + i^{M) =a (36) 

From this, we get a Hamiltonian for the variable 1Z: 



w<2) = Jm p «' < 37 > 



where P-r is a conjugate momentum of 1Z. The equation of the motion for 1Z becomes 

K aa + ( 3 + ^ + 2 A In ft Q = 0. (38) 



H da \H 

By introducing a new variable Q = 4> a 1Z = <p ao 5a + (fic$C, Q satisfies 
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Qm + (3 + Q„ + _ _L_ (e^tf^) = 0. (39) 

This is nothing but a Mukhanov equation for the long- wavelength limit. For 0-mode pertur- 
bation, the only gauge freedom is infinitesimal time coordinate transformation t — > t + St. 
Then variables transform as Sa — > Sa + aSt, S(p — > S(p + (j)St. A possible gauge invariant 
combination of these variables is Sa — (a/4>)S<f> and this has a meaning of perturbation of 3 
curvature in co-moving gauge (intrinsic curvature perturbation). In our calculation, we use 
a as a time parameter and this is equivalent to using co-moving gauge. Therefore the vari- 
able 71 is the gauge invariant intrinsic curvature perturbation or the Bardeen's parameter 
in the long- wavelength limit. 

The two independent solutions are given by 

Q = K,Ac- (40) 

The explicit form of these solutions can be obtained using Eq.(27), (28), (29). Assuming 
that the lapse iV is a function of the scalar field 0, we have 

1 ft d(t> 

t+t " = -2l inr; (41) 

Regarding as a function of (t + t , C) and differentiating the both side of this equation 
with respect to t and C, 



d 



1 



dt 2NH<p 
d 4>c [* H 4>c 



bc- = -wh: + J wm^- (42) 



Choosing the scale factor as a time parameter(t = a, N = 1/H), 

Ha, 



<Pa = - 2 



H r " ; 

r<t> H^n r a e^ 3(a+Qo) 

NH *J iv# rf0 = " 60Q i ~tfjT da ' (43) 



We therefore obtained the gauge invariant 0-mode solution by differentiating the background 
solution with respect to two parameters a and C. 
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B. two scalar fields case 



We consider a flat FRW universe with two scalar fields 0i and 2 . The Hamilton- Jacobi 
equation is 

- § + + + e 6a V(<j> u 2 ) = 0. (44) 

The background solution is obtained by = — 2e 3o if(0i, 2 ), 

- 3H 2 + 2Hl + 2Hl + Vfa, 2 ) = 0, (45) 
£ = H, (46) 
| = -2H^, | = -2^ 2 - (47) 

The solution of Eq.(45) has two constants of integration Ci,C 2 and can be written H = 
H((f>i,<f>2,Ci,C2)- By differentiating the H-J equation (45) with respect to these constants 
and using the evolution equation (46) and (47), it can be shown 

e 3a H cl = D 1 , e 3a H C2 = D 2 , (48) 

where D 1: D 2 are constants. Let us define new constants a and / by D ± — e~ 3a °,D 2 = 
e~ 3a "f. Then 

H C i = e- 3{a+ao \ (49) 

Hc2 



Hci 



f- (50) 



From this we have 



a + a — 2 , Ci, C 2 ), (51) 
f = f(<f> 1 ,cf> 2 ,C 1 ,C 2 ). (52) 

ao is the origin of scale factor. If C±, C 2 are fixed, / determines trajectories in configuration 
space (0x, 2 ). We split the generating function into the background and the perturbed part: 
S = + F. F satisfies 
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S * ] P _i_ P _i_ P i F <h 

Change variables from (0 1? 2 , &) to (a, a , C±, f, C 2 ): 

d 



0. 



d_ 

da 
d 

Wi 

d 



dao 



r 9 d 



<9CV 



dd 

~ f d d \ ^ d . d „ d 



(9a <9a 



d 

-r + 



d 



902 \ da dao , 

The H-J equation becomes 



+ C U 



d 



df 
d 



dC 2 
d 



dCl + Uz d f +C ^ d C 2 



(53) 



(54) 
(55) 
(56) 



da 12 V da 



dF dF dF 



<9C 2 



1/ /9F dF\ dF dF dF' 



1 / /<9F <9F\ 



<9F dF dF" 



To extract a linear perturbation part, we replace ^ 

9 



9 9 



95/ ' dCi 



(57) 



9 9 



9<5C*i ' 9C 2 



9ao 9<5ao ' 9/ 

a5C , 2 . Using a new variable 71 = 5ao + aciSCi + <Sc2<5C2 and assuming a quadratic form of 
F with respect to 7£ and Sf, we have 



dF 



dF 



e3aH: da + 2{ a ^dn + f ^d6f 



dF 



dF" 
d6f, 



= 0. 



(58) 



Define 2x2 matrices 



X 



Via, <Plf 
02a, 02/ 



The H-J equation becomes 



<9F 1 

&5 + 2 



9£ dF_ 

97?. 95/ 



9£ 
d1l 

dF 
dSf 



0. 



We get a Hamiltonian for the perturbation variables (1Z, Sf): 



(59) 



(60) 



W (2) 



PR 
^5/ 



(61) 
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where Pn and P$f are conjugate momentum of 1Z and 5f, respectively. The evolution 
equation becomes 



K 



+ A (A- 1 



■R 



o, 



(62) 



Introducing a new variable 



- aa 










Qi 




K 


Q = 


= X 






Q 2 







(63) 



Q satisfies 



Qnn (3+^) Q n + UQ = (). 



(64) 



where M is a 2 x 2 matrix defined by 



(65) 



This is the Mukhanov equation for two scalar field system in the long- wavelength limit [8]. 
The four independent solutions of this equation are 



Q = 


01 


) 


01 




01 




01 




02 




02 


/ 


02 


Ci 


02 



(66) 



c 2 



These solutions can be written in another form: 



Q = X 



d 3 

G?4 



(67) 



where ^1,^2,^3,^4 are arbitrary constants. Using the Mukhanov's variable, 1Z and Sf can 
be expressed as 



n 






= X~ X Q = 


Sf 





U1Q1 + ffoQz 

In two scalar field case, using a as a time parameter means 



(68) 
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5a = a^S^i + a</, 2 502 = 0. (69) 

From the relation 

dadH H 

a * = Md^ = ~^rJ^ (70) 

we have 5a = J2i 4>i^4>i = 0- So using a as a time parameter is equivalent to using co- 
moving gauge and the variable 1Z is the gauge invariant intrinsic curvature perturbation in 
the long- wavelength limit. 

We can rewrite the expression of 7Z using Eq.(70): 

n=- 1 ^-{<p la Q l + <p 2a Q 2 ). (71) 

This gives a relation between the curvature perturbation 1Z which corresponds to the 
Bardeen's parameter ( and the Mukhanov's variable Qi,Q 2 in the long- wavelength limit. 
Qi,Q 2 are given by Eq.(67) which is written by using background solution. The important 
point is that we have not assumed any approximation to derive the long-wavelength solution. 
We can get the gauge invariant gravitational potential <£> by considering perturbation of the 
evolution equation (46): 

5a - H5N = H^Sfa. (72) 

i 

In longitudinal gauge, 5N = —5a = $ and Qi = 5<pi + (4>i/H)^. The above equation 
becomes 

$ + H$ = \Y,j>Mi- ( 73 ) 

The curvature perturbation 71 can be written using $: 

Iti j 

TT 

= $--^($ + HQ). (74) 



Therefore we have 



$ = "I / ^IP 11 (75) 
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C. application to inflationary model 



We demonstrate our 0-mode solution (66) of perturbation reproduces previously obtained 
results [10-12] in inflationary model using the slow rolling approximation. If we concentrate 
on growing mode solutions, the gauge invariant 0-mode solution is given by 



(76) 



where A, B are arbitrary constants. Substituting this solution to the expression of the 
curvature perturbation (71), we have 



Qi 


= A 


01a 


+ B 


0i/ 






Q 2 




02a 




02/ 



n = a-b 



H 

2H« 



Ha01/ + 02a02/) 



(77) 



The first term is the contribution of adiabatic fluctuations and remains constant in time. 
The second term is the contribution due to iso-curvature or entropic fluctuations which 
arises whenever there are two or more scalar fields. Using the slow rolling approximation, 
H-J equation and evolution equations become 



iJ 2 ^i\/(0 1 ,0 2 ), 

2„ _ 
~H H * X ~ ' 







V, 



la 



01 



02a = 



-H, 



V ' 

v„ 



(78) 



H -4>2 y 

In this approximation, we do not have parameters C\,Ci contained in the Hubble function 
H and we neglect decaying mode solutions of perturbations. In the inflationary phase, 
H/H <C 1 and the gauge invariant gravitational potential becomes 

H 



4>2 



2F 2 



(79) 



For general form of the potential \^(0i,0 2 ), we can not write down the solution of the 
evolution equation (78). We consider two examples of the potential in which case we have 
integral form of the slow rolling solution. 

Case 1 : V((j) 1: 2 ) = V"i(0i) + V 2 (02)- By integrating the evolution equation (78), we have 
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/V f V 

d(j)i—A / d&jr^-, (80) 

J V l,4>l •> V 2,<f>2 

where ao, f are constants of integration. By differentiating the above solution with respect 
to the parameter /, we have 

*» = -wh v ^- (82) 

Substituting these solutions to the expression of 1Z (77), we have 

The gauge invariant gravitational potential is 

H a BV la V 2 -V 2a V 1 , , 

$ = A— ±£i_i. (84) 



Case 2 : l / (0i,0 2 ) = Vi(0i)V2(02)- By integrating the evolution equation (78), we have 

a + a = - f d^-^- f i 



4>i •> 2V2,02 ' 



/ = / # * - I d<p", (85) 



/ WITT- 

where ao, / are constants of integration. By differentiating the above equations with respect 
to /, we have 

, _ Vi,<f>i 
<h,f ~ ^ 

Substituting these solutions to the expression of 1Z (77), we have 



(86) 



TZ = A-- [vi) 2 [v > } 2 . (87) 

2 m+m 



The gauge invariant gravitational potential is 
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IV. SUMMARY 



We have shown that the solution and the evolution equation of the gauge invariant pertur- 
bation in the long-wavelength limit can be derived within the mini-super space model. The 
solution is obtained by differentiating the background solution with respect to parameters 
contained in the background system. These parameters completely specify the background 
solutions. If we need to know the behavior of the long-wavelength perturbation only, this 
method is very useful because we does not have to solve the perturbation equation. What 
we have to do is to solve the background system and identify the parameters contained in 
the background solution. The analysis of the long-wavelength perturbation in the model 
that involves multiple scalar field becomes easier. 

The most important feature of our 0-mode solution is it does not rely on the slow roll 
approximation which is usually assumed to analyze the perturbation in the multi-field infla- 
tionary model [10-13]. As an application, a model of reheating after inflation was considered 
[9]. The system consists of a oscillating inflaton field and a massless scalar field which inter- 
acts with the inflaton. By using the approximation which can include non-linear effect, we 
obtained background solution and the behavior of the long-wavelength curvature perturba- 
tion was obtained using the method presented in this paper. 

As a straightforward extension, it is possible to treat spatially closed and open FRW 
universe and anisotropic Kasner type universe. Instead of scalar fields, perfect fluid can be 
included. These subjects are left as a future exercise. 

This work is partly supported by the Grant-In-Aid for Scientific Research of the Ministry 
of Education, Science, Sports and Culture of Japan(09740196). 
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